Numerical integration of differential equations in Matlab